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Abstract 



Let (l)Rh = f, < x < L, Rh = J Q L R(x,y)h(y) dy, where the kernel R(x,y) 
satisfies the equation QR = P5(x — y). Here Q and P are formal differential oper- 
ators of order n and m < n, respectively, n and m are nonnegative even integers, 
n > 0, m > 0, Qu := q n (x)u^ + YTjZl qj(x)u®, Ph := + Pj(x)hV> , 

Qn(x) > c > 0, the coefficients qj(x) and Pj(x) are smooth functions defined on R, 
5(x) is the delta-function, / £ H a (0,L), a := H a is the Sobolev space. 

An algorithm for finding analytically the unique solution h £ H~ a (0,L) to 
(1) of minimal order of singularity is given. Here H~ a (0,L) is the dual space to 
H a (0, L) with respect to the inner product of L 2 (0, L). 

Under suitable assumptions it is proved that R : H~ a (0,L) — > H a (0,L) is an 
isomorphism. 

Equation (1) is the basic equation of random processes estimation theory. Some 
of the results are generalized to the case of multidimensional equation (1), in which 
case this is the basic equation of random fields estimation theory. 



1 Introduction 



In monograph pQ (see also [3]) estimation theory for random fields and processes is 
constructed. The estimation problem for a random process is as follows. Let u(x) = 
s(x) + n(x) be a random process observed on the interval (0, L), s(x) is a useful signal 
and n(x) is noise. Without loss of generality we assume that s(x) = n(x) = 0, where the 
overbar stands for the mean value, u*(x)u(y) := R(x,y), R(x,y) = R(y,x), u*(x)s(y) : = 
f(x, y), and the star here stands for complex conjugate. The covariance functions R(x, y) 

Key words: integral equations, estimation theory 

Math subject classification: 45A05, 45H05, 93E10, 93E11, 34A30 



and f(x, y) are assumed known. One wants to estimate s(x) optimally in the sense of 
minimum of the variance of the estimation error. More precisely, one seeks a linear 
estimate 

Tu= h(x,y)u(y)dy, (1.1) 
Jo 

such that 



\(Tu)(x) — s(x)\ 2 = min. (1.2) 

This is a filtering problem. Similarly one can formulate the problem of optimal estimation 
of (As)(x), where A is a known operator acting on s(x). If A = I, where I is the identity 
operator, then one has the filtering problem, if A is the differentiation operator, then one 
has the problem of optimal estimation of the derivative of s, if As = s(x + Xq), then one 
has an extrapolation problem, etc. The kernel h(x, y) is, in general, a distribution. As 
in p, one derives a necessary condition for h to satisfy (1.2): 

L 

R(x,y)h(y,z)dy = f(x,z), 0<x,z<L. (1.3) 
Since z enters as a parameter in (1.3), the basic equation of estimation theory is: 

Rh := [ R(x,y)h(y)dy = f(x), < x < L. (1.4) 



The operator in L 2 (0,L) defined by (1.4) is symmetric. In |1] it is assumed that the 
kernel 

R(x, y) = J™ V, A) dp(X), (1.5) 

where P(X) and Q(X) are positive polynomials, $>(x,y,X) and dp(X) are spectral kernel 
and, respectively, spectral measure of a selfadjoint ordinary differential operator £ in 
L 2 (M), degQ(A) = q, degP(A) = p < n, p > 0, ordi := a > 0. Actually in 1J the 
multidimensional case, when x, y G M r , r > 1, and £ is a selfadjoint elliptic operator in 
L 2 (W) is also considered. 

It is proved in |T] that the operator R : H~ a (0,L) — > H a (0,L), a := "^^a is an 
isomorphism. By H a (0,L) the Sobolev space W a,2 (0,L) is denoted, and H~ a (0,L) is 
the dual space to H a (0, L) with respect to L 2 (0, L) := H°(0, L) inner product. Namely, 
H~ a (0, L) is the space of distributions h which are linear bounded functionals on H a (0, L). 
The norm of h G H~ a (0, L) is given by the formula 

lls.ll l^'^l (i a\ 

\\ h \\H-«(0,L) = SU P Fl ' ( L6 ) 

g£H<x(0,L) \\9\\H a {Q,L) 

where (h, g) is the L 2 (0, L) inner product if h and g belong to L 2 (0, L). 

One can also define H~ a (0, L) as the subset of the elements of H~ a (M>) with support 
in [0,L]. 
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In this paper we generalize the class of kernels R(x,y) defined in (1.5): we do not 

use the spectral theory, do not assume t to be selfadjoint, and do not assume that the 
operators Q and P commute. 

We assume that 

QR = P5(x-y), (1.7) 

where Q and P are formal differential operators of orders n and m respectively, n > m > 
0, n and m are even integers, 8(x) is the delta-function, 

n m—1 

Qu:=^2qj(x)u^, q n {x) > c> 0, Ph := /i (m) + ^ Pj(x)h®, (1.8) 

3=0 3=0 

where qj and pj are smooth functions defined on R. We also assume that the equation 
Qit = has | linearly independent solutions uj G L 2 (— oo, 0) and | linearly independent 
solutions G L 2 (0,oo). In particular, this implies that if Qh — 0, h G i/ a (R), a > 0, 
then /i = 0, and the same conclusion holds for h G iJ^R) for any fixed real number (3, 
including negative f3, because any solution to the equation Qh = is smooth: it is a 
linear combination of n linearly independent solution to this equation, each of which is 
smooth and none belongs to L 2 (R). 

Let us assume that R(x, y) is a selfadjoint kernel such that 

CilMI 2 . < (R<p, <p) < C2IMI 2 ., ci = const > 0, Vc/? G C °°(R), (1.9) 

where (•, •) is the L 2 (R) inner product, \\(p\\- '■— ll^llff-"^) := a '■= IMI/3 := 

1 1 y 1 1 -ff" (M) j an( i we use below the notation \\(p\\+ '■= \\<p\\H a (o,L) '■= \\^p\\h+- The spaces 
H a (0, L) and H~ a (0, L) are dual of each other with respect to the L 2 (0, L) inner product, 
as was mentioned above. If (p G H~ a (0,L), then y> G iif _a (R), and the inequality (1.9) 
holds for such (p. By this reason we also use (for example, in the proof of Theorem 1.1 
in Section 2) the notation H~ for the space H~ a (0, L). 

Assumption (1.9) holds, for example, if the following inequalities (1.10) and (1.11) 
hold: 

C3|M|-a+n < HQVH-a < C 4 1| if\\ - a +n, C 3 , C 4 = Const > 0, G C~(M), (1-10) 

cslMlL™ < (PQ*<p,<p) <c 6 ||^|| 2 1±22 , V^gC °°W, (1.11) 
2 2 

where is a formally adjoint to Q differential expression, and C5 and Cq are positive 

constants independent of ip G C^°(R). The right inequality (1.11) is obvious because 

ord PQ* = n + m, and the right inequality (1.10) is obvious because ordQ* = n. 

Let us formulate basic results of this paper. 

Theorem 1.1. If (1.9) holds, then the operator R, defined in (1.5), is an isomorphism 
ofH- a (0,L) onto H a (0,L), a= 
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Theorem 1.2. // (1.7), (1.10) and (1.11) hold, then (1.9) holds and R : H~- a (0,L) -> 
H a (0,L) is an isomorphism. 

Theorem 1.3. If (1.7), (1.10) and (1.11) hold, and f G H n (0,L), then the solution to 
(1-4) in H~ a (0, L) does exist, is unique, and can be calculated analytically by the following 
formula: 

px n—ot—1 

h= G(x,y)Qfdy + V [aj(—l) j G^\x,0) + af(—l) j G^\x, L)] , (1.12) 
Jo j=0 

where are some constants and G(x, y) is the unique solution to the problem 

PG = 6{x-y), G(x,y) =0 for x <y. (1.13) 

The constants are uniquely determined from the condition h(x) = for x > L. 

Remark 1.4. The solution h G H~ a (0,L) is the solution to equation (1.4) of minimal 
order of singularity (see ff$). 

Remark 1.5. If P = 1 in (1.7) then the solution h to (1-4) of minimal order of singular- 
ity, h G i/~2 (0, L), can be calculated by the formula h = QF, where F is given by (2.9) 
(see below) and u + and U- are the unique solutions of the problems Qu + = if x > L, 
u [ i\L) = f {j) (L), < j < § - I, u+ (oc) = 0, and Qu„ = if x < 0, u U) (0) = f {j) (0), 
<j < f -1, u_(-oo) = 0. 



2 Proofs 

Proof of Theorem 1.1. The set C^°(0, L) is dense in H~ a (0, L) (in the norm of H~ a 
Using the right inequality (1.9), one gets: 



\R\\h-->h+ = sup ' < ca, (2.1) 
heH- \\' l \\- 



{Rh, h) 

\n\\fj-^H+ — sup 
heH~ 

by the symmetry of R in L 2 (0, L) (cf [TJ p. 40]). This implies !#---*.#+ < ci- Using the 
left inequality (1.9), one gets: ci||/i|| 2 < + so 

ci||ft||_ < \\Rh\\ + . (2.2) 

Therefore 

WR-'Wh+^h- < -■ (2.3) 
c% 

Consequently, the range Ran(R) of R is a closed subspace of H + . In fact, Ran(R) = H + . 

Indeed, if Ran(R) ^ H + , then there exists a g G H~ such that = (Ri/j,g) Vip G H~ . 

Taking i/j = g and using the left inequality (1.9) one gets \\g\\- = 0, so g = 0. Thus 
Ran(R) = H+. 

Theorem 1.1 is proved. □ 
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Proof of Theorem 1.2. From (1.7) and (1.8) it follows that the kernel R(x,y) defines a 
pseudodifferential operator of order —2a = m — n. In particular, this implies the right 
inequality (1.9). In this argument inequalities (1.10) and (1.11) were not used. 

Let us prove that (1.10) and (1.11) imply the left inequality (1.9). 

One has 

HQVIU < C|MU_ Q) W>gC °°(R), (2.4) 

because ordQ* = n. Inequality (1.10) reads: 

c 3 |MI-a+n < ||QV||_ Q < c 4 ||<p||_ a+n , VVeC£°(R), (2.5) 

where c 3 and c 4 are positive constants. If (2.5) holds, then Q* : H~ a+n (R) — > H~ a (R) is 
an isomorphism of H~ a+n (Ei) onto H~ a (R) provided that N(Q) := {w : Qw = 0, w G 
if Q (R)} = {0}. Indeed, if the range of Q* is not all of H~ a (R), then there exists an 
w ^ 0, w G H a (R) such that (Q*<p,w) = G C °°(M), so Qw = 0. If Qw = and 
w E H a (R), then, as was mentioned below formula (1.8), it follows that w — 0. This 
proves that Ran(Q*) = #~ Q (M). 

Inequality (1.11) is necessary for the left inequality (1.9) to hold. Indeed, let ip = Q*tp, 
up E C °°(R), then (1.9) implies 

c 5 |MI 2 _ a+n < c\\Q*<p\\ 2 _ a < (i2Q>, QV) = (QRQ*^ <p) = (PQV, ^), (2.6) 

where c > here (and elsewhere in this paper) stands for various estimation constants. 
Because — a + n = ^±^ ; inequality (2.6) is the left inequality (1.11). The right inequality 
(1.11) is obvious because the order of the operator PQ* equals to n + m. 

Let us prove now that inequalities (1.11) and (1.10) are sufficient for the left inequality 
(1.9) to hold. 

Using the right inequality (1.10) and the left inequality (1.11), one gets: 
c||^|| 2 _ a <c 5 || V 9||i ±2ri <(PQV^) = W,V0, V> = QV, V^eC °°(R). (2.7) 

2 

Let us prove that the set {ip} = {QVlv^ec^OR) is dense in H~ a (0,L). Assume the 
contrary. Then there is an h G H~ a (0, L), h ^ 0, such that (Q*ip,h) = for all 
if G Cg°(R). Thus, (<p,Qh) = for all tp G C^°(R). Therefore Q/i = 0, and, by the 
argument given below formula (1.8), it follows that h — 0. This contradiction proves 
that the set {Q*^}\/ ipeC g c (R) is dense in H~ a (0,L). 

Consequently, (2.7) implies the left inequality (1.9). The right inequality (1.9) is an 
immediate consequence of the observation we made earlier: (1.7) and (1.8) imply that R 
is a pseudodifferential operator of order —2a = —(n + m). 

Theorem 1.2 is proved. □ 

Proof of Theorem 1.3. Equations (1.4) and (1.7) imply 

Ph = g:= QF. (2.8) 
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Here 

x < 0, 

/, < x < L, (2.9) 
u + , x > L, 

where 

Qw_ = 0, x<0, (2.10) 

Qu + = 0, x > L, (2.11) 
and U- and u + are chosen so that F e iJ a (IR). This choice is equivalent to the conditions: 

«i?' ) (0) = /W(0) J 0<j<a-l, (2.12) 

U J' ) (L) = /^(L), 0<j<a-l. (2.13) 
If F G # a (M), then g := QF G if a_n (R) = and, by (2.9), one gets: 

n— a— 1 

£ = Q/+ £ [a-^(x) + a+^(x-L)], (2.14) 
j=o 

where af are some constants. There are n — a = !1 ^ n - constants a+ and the same number 
of constants aj. 

Let G(x, y) be the fundamental solution of the equation 

PG = 6(x - y) in R, (2.15) 

which vanishes for x < y: 

G(x,y) = for x < y. (2.16) 
Claim. Such G(x, y) exists and is unique. It solves the following Cauchy problem: 



PG = 0, x>y, G^(x,y) 
satisfies condition (2.16), and can be written as 



= Sj, m -i, 0<j<m-l, (2.17) 

x=y+0 



G(x,y) = ^2cj{y)(pj{x), x>y, (2.18) 

where fj(x), 1 < j < m, is a linearly independent system of solutions to the equation: 

Ptp = 0. (2.19) 
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Proof of the claim. The coefficients Cj(y) are defined by conditions (2.17): 



3=1 



m—l-i 



< k < m- 1. 



(2.20) 



The determinant of linear system (2.20) is the Wronskian W(tpi, . . . , <p m ) ^ 0, so that 
Cj(y) are uniquely determined from (2.20). 

The fact that the solution to (2.17), which satisfies (2.16), equals to the solution to 

(2.15) - (2.16) follows from the uniqueness of the solution to (2.15) - (2.16) and (2.17) 
- (2.16), and from the observation that the solution to (2.15) - (2.16) solves (2.17) - 

(2.16) . The uniqueness of the solution to (2.17) - (2.16) is a well-known result. 

Let us prove uniqueness of the solution to (2.15) - (2.16). If there were two solutions, 
Gi and G 2 , to (2.15) - (2.16), then their difference G := G\ — G 2 , would solve the 
problem: 

PG = in R, G = for x < y. (2.21) 

By the uniqueness of the solution to the Cauchy problem, it follows that G = 0. Note 
that this conclusion holds in the space of distributions as well, because equation (2.21) 
has only the classical solutions, as follows from the ellipticity of P. Thus the claim is 
proved. □ 

From (2.8) and (2.14) - (2.16) one gets: 



(*X (*X Tl ft 1 

h= G(x,y)Qfdy + G(x,y) V [aj 6® (y) + a+ 6® (y - L)] dy 

JO JO „_ n 

G(x,y)Qfdy 



j=o 



n—a—l 



f 

JO 



E (-v 

3=0 



G^(x,y) aJ + G^(x,y) 
y=o 



y=L 



G(x,y)Qfdy + H(x). 



It follows from (2.22) that h 6 H- a (R) and 

h = for x < 0, 



^2.22) 



(2.23) 



that is, (h, (f) = Vy? G C™(M) such that supp ip C (— oo, 0). In order to guarantee that 
h E H~ a (0, L) one has to satisfy the condition 



h = for x > L. 



(2.24) 



Conditions (2.23) and (2.24) together are equivalent to supp h C [0, L\. Note that al- 
though Qf E i/ _2L ^ 2i (0, L), so that Qf is a distribution, the integral J* G(x, y)Qf dy = 
X^oo G(x, y)Qf dy is well defined as the unique solution to the problem Pw = Qf, w = 
for x < 0. 
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Let us prove that conditions (2.23) and (2.24) determine the constants af, < j < 
— 1, uniquely. 

If this is proved, then Theorem 1.3 is proved, and formula (2.22) gives an analytical 
solution to equation (1.4) in H~ a (0, L) provided that an algorithm for finding af is given. 
Indeed, an algorithm for finding G(x, y) consists of solving (2.16) - (2.17). Solving (2.16) 
- (2.17) is accomplished analytically by solving the linear algebraic system (2.20) and 
then using formula (2.18). We assume that m linearly independent solutions tpj(x) to 
(2.19) are known. 

Let us derive an algorithm for calculation of the constants af, < j < — 1, from 
conditions (2.23) - (2.24). 

Because of (2.16), condition (2.23) is satisfied automatically by h defined in (2.22). 
To satisfy (2.24) it is necessary and sufficient to have 

L 

G(x,y)Qfdy + H(x)=Oiorx>L. (2.25) 

By (2.18), and because the system {<£j}i<j< m is linearly independent, equation (2.25) is 
equivalent to the following set of equations: 

n-\-m 

/ c k (y)Qfdy+ 

Let us check that there are exactly m independent constants af and that all the constants 
af are uniquely determined by linear system (2.26). 

If there are m independent constants af and other constants can be linearly repre- 
sented through these, then linear algebraic system (2.26) is uniquely solvable for these 
constants provided that the corresponding homogeneous system has only the trivial so- 
lution. If / = 0, then h — 0, as follows from Theorem 1.1, and g = in (2.14). Therefore 
af = Vj, and system (2.26) determines the constants af Vj uniquely. 

Finally, let us prove that there are exactly m independent constants af. Indeed, in 
formula (2.8) there are | linearly independent solutions uj G L 2 (— oo,0), so 

n 

u -=i2 b i u i' ( 2 - 27 ) 

and, similarly, u + in (2.8) is of the form 

n_ 

2 

u + = Y,b J uf, (2.28) 
where uf G L 2 (0, oo). Condition F G H a (R) implies 

n 
2 

b-(uj)W = at x = 0, 0<A;<a-l = - 1, (2.29) 



Jo) 



(OK 



= 0, 1 < k < m. 



(2.26) 
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and 

n 

^ + «) (fc) = f {k) at x = L, < k < - 1. (2.30) 

3=1 

Equations (2.29) and (2.30) imply that there are | — = y independent constants 6~ 
and — independent constants and the remaining n — m constants b~ and b1~ can be 
represented through these m constants by solving linear systems (2.29) and (2.30) with 
respect to, say, first !i y 2; constants, for example, for system (2.29), for the constants bj , 
1 < j < ^^y 1 . This can be done uniquely because the matrices of the linear systems 
(2.29) and (2.30) are nonsingular: they are Wronskians of linearly independent solutions 
{^j }i<j< n ~ m and {uj }i<j< "- m - 

The constants can be expressed in terms of b^ and / by linear relations. Thus, 
there are exactly m independent constants aj. This completes the proof of Theorem 1.3. 
□ 

Remark 2.1. In JU pp. 317-337] a theory of singular perturbations for the equations of 
the form 

eh £ + Rh £ = f (2.31) 

is developed for a class of integral operators with a convolution kernels R(x, y) = R(x — y). 
This theory can be generalized to the class of kernels R(x, y) studied in the present paper. 
The basic interesting problem is: for any e > equation (2.31) has a unique solution 
h £ G L 2 (0,L); how can one find the asymptotic behavior of h e as e —>■ 0? The limit 
h of h £ as e — > should solve equation Rh = f and, in general, h is a distribution, 
h G H~ a (0, L) . The theory presented in pp. 317-337] allows one to solve the above 
problem for the class of kernels studied in the present paper. 

Remark 2.2. Theorems 1.1 and 1.2 and their proofs remain valid in the case when 
equation (1-4) is replaced by the 

Rh:= / R(x,y)h{y)dy= f, x G D. (2.32) 
Jd 

Here D C W , r > 1, is a bounded domain with a smooth boundary S, D is the closure 
of D, R(x,y) solves (1.7), where P and Q are uniformly elliptic differential operators 
with smooth coefficients, ordP = m > 0, ordQ = n > m, equation Qh = has only the 
trivial solution in H^(W) for any fixed real number j3. Under the above assumptions, one 
can prove that the operator defined by the kernel R(x, y) is a pseudodifferential elliptic 
operator of order —2a, where a := We do not assume that P and/or Q are 

self adjoint or that P and Q commute. Therefore Theorems 1.1 and 1.2 for equation (2.32) 
generalize the results from *!]. An analog of Remark 2.1 holds for the multidimensional 
equation (2.31) as well. Equation (2.32) is the basic integral equation of random fields 
estimation theory £7]/. 
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In the short announcement by L.Piterbarg, Differentsialnye Uravneniya, 17, (1981), 
2278-2279, (in Russian), some results concerning equation (2.32) are announced, but the 
author is not aware of the published proofs of these results. 
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